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Energy dissipation in magnetohydrodynamic (MHD) turbulence is known to be highly inter¬ 
mittent in space, being concentrated in sheet-like coherent structures. Much less is known about 
intermittency in time, another fundamental aspect of turbulence which has great importance 
for observations of solar flares and other space/astrophysical phenomena. In this Letter, we 
investigate the temporal intermittency of energy dissipation in numerical simulations of MHD 
turbulence. We consider four-dimensional spatiotemporal structures, “flare events”, responsible 
for a large fraction of the energy dissipation. We find that although the flare events are often 
highly complex, they exhibit robust power-law distributions and scaling relations. We find that the 
probability distribution of dissipated energy has a power law index close to a « 1.75, similar to 
observations of solar flares, indicating that intense dissipative events dominate the heating of the sys¬ 
tem. We also discuss the temporal asymmetry of flare events as a signature of the turbulent cascade. 


PACS numbers: 52.35.Ra, 95.30.Qd, 96.60.lv, 52.30.Cv 


Introduction .— Intermittency plays a major role in tur¬ 
bulence by causing processes such as energy dissipation 
and particle acceleration to be highly localized in coher¬ 
ent structures. It also forestalls efforts toward a complete 
theory of turbulence. Many tools have been employed to 
study intermittency, including structure functions |1, 2], 
scale-dependent kurtosis Q, topological methods [j. 5|, 
and statistics of discontinuities [61, 3]. However, past 
studies have mainly focused on spatial intermittency, giv¬ 
ing limited information about the dynamics. In order 
to understand the temporal aspects of intermittency, in¬ 
cluding characteristic timescales of structures as well as 
their interactions and stability, a broader framework is 
needed. 


A promising new paradigm is the statistical analysis 
of coherent structures, which is robust and informative 
for studies of intermittency. The occurrence rates, in¬ 
tensities, and morphology of structures yield insight to 
the inhomogeneity, anisotropy, and characteristic scales 
of the dynamics. Coherent structures can be simply iden¬ 
tified as regions in space bounded by an isosurface of 
some field. This was used to study vorticity filaments 
in hydrodynamic turbulence magnetic structures 

in the kinematic dynamo [ll| , and dissipative structures 
in magnetohydrodynamic (MHD) turbulence 12 -0] and 
ambipolar diffusion MHD |l6[. Since coherent structures 
and intense dissipative events are experimentally observ¬ 
able, there are many practical applications including so¬ 
lar flares, instabilities in fusion devices Il7j . and radiative 
signatures in optically thin astrophysical plasmas, e.g., in 
black-hole accretion disk coronae [l8|, hot accretion flows 
0, and jets [0; in pulsar wind nebulae [e.g. mm ; 
and possibly in the hot gas in galaxy clusters. 

This Letter addresses some fundamental aspects of in¬ 


termittency in MHD turbulence. A major question is 
whether, in the limit of large Reynolds number, energy 
dissipation is dominated by a few intense, large-scale 
events or by many weak, small-scale events. A related 
question is whether there is an inherent relationship be¬ 
tween spatial intermittency and temporal intermittency, 
e.g., whether larger structures retain their coherency in 
time. These temporal aspects of intermittency have been 
practically unexplored in previous MHD studies. 

In this Letter, we extend a framework previously de¬ 
veloped for the statistical analysis of dissipative struc¬ 
tures 0 into the temporal realm, thereby considering 
4D spatiotemporal objects representing flare events. We 
apply this novel methodology to study intermittency in 
numerical simulations of strong incompressible MHD tur¬ 
bulence. We describe the distributions, scalings, and evo¬ 
lution of flare events by characterizing their length scales, 
durations, dissipated energies, and peak energy dissipa¬ 
tion rates. These are the first results on the fundamental 
properties of the combined spatial and temporal inter¬ 
mittency of energy dissipation in 3D MHD turbulence. 


The primary questions addressed here for MHD turbu¬ 
lence are also fundamental for the solar corona. In fact, 
our approach has strong similarities with observational 
studies of solar flares 23- 32| and stellar flares 33-0, 
which use the time-series of X-ray and extreme UV emis¬ 
sions to measure the duration, peak intensity, and fluence 
of flares, from which dissipated energy is inferred. For 
the solar corona, a measurement of central importance 
is the probability distribution for dissipated energy, due 
to its role in assessing the nanoflare model for coronal 
heating 37, 38j. This distribution exhibits a power law 
over eight orders of magnitude, with an index near —1.8, 
somewhat shallower than the critical index of —2 required 
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for nanoflares to dominate the overall heating [39|. 

We compare our results with the observed statistical 
properties of solar flares. We stress that there are several 
basic differences between incompressible MHD turbu¬ 
lence and the solar corona. In contrast to volumetrically- 
driven turbulence, the solar corona is modeled by force- 
free MHD with slowly-driven, line-tied boundaries. Fur¬ 
thermore, kinetic effects may become important during 
reconnection. Despite these differences between our sim¬ 
ulations and the solar corona, we find that the statisti¬ 
cal properties of flare events have multiple similarities in 
both cases. This suggests that MHD turbulence may play 
a role in the energetics of the corona [30, Eo| , a possibil¬ 
ity that should be investigated more carefully in future 
studies. 


Method .— We perform simulations of reduced MHD, 
applicable since the uniform background magnetic field 
Bo = Bqz is strong relative to turbulent fluctuations, 
B 0 « 56 rms . The equations are [4l!j 

(Jf t TV A - V,| + z* • Vx) z± = -V±P + + ff 

V_L-Z ± =0 (1) 

where z ± = v ± b are the Elsasser variables (perpen¬ 
dicular to Bo), v is the velocity field, b is the fluc¬ 
tuating magnetic field (in units of the Alfven velocity, 
Va = Bq/^/Attpq, where po is the uniform plasma den¬ 
sity), P is the total pressure, and /* is the external 
forcing. We use uniform fluid viscosity v and magnetic 
diffusivity rj with v = rj. We consider structures in the 
current density j = j z = z ■ Vi x 6; the resistive energy 
dissipation rate per unit volume is rjj 2 . 

Equations © are solved using a fully dealiased 3D 
pseudo-spectral algorithm (for details, see J42|). The pe¬ 
riodic box is elongated in z by a factor of L»/Lj_ = 6, 
where L± = 2n is the box size in simulation units. 
Timescales are in units of eddy turnover times, r e ddy = 

L x/(27ru rms ) ~ 1. The turbulence is driven at large 
scales by colliding Alfven waves, generated from statis¬ 
tically independent random forces at wave-numbers 
27 t/L± < k XtV < 2(27 t/L±), k z = 2tt/L\\. The forc¬ 
ing is solenoidal, has random Fourier coefficients taking 
Gaussian values that are refreshed independently approx¬ 
imately 10 times per eddy turnover time, and has ampli¬ 
tude such that 6 rms ~ u rms ~ 1. The Reynolds number 
is given by Re = v Tms (L±/ 2-k)/v. 

We analyze snapshots dumped at a cadence (At) -1 
from four simulations shown in Table |U The main results 
are from runs with 512 3 resolution, with the Re = 1250 
case having the highest cadence and longest time interval. 
Due to computational constraints, At is larger than the 
internal time step in the simulation. The minimum ca¬ 
dence required to properly track structures is estimated 
by requiring that the distance advected by the flow dur¬ 
ing At is less than the typical current sheet thickness, 


giving u rms A t < & r ms/jthr- The cadences are compara¬ 
ble to this value and the results show convergence with 
cadence. 


TABLE I. Simulations 


Sim. 

Res. 

Re 

At- 1 

Time interval 

1 

256 3 

800 

64 

10.0 

2 

512 3 

800 

32 

12.2 

3 

512 3 

1250 

64 

15.6 

4 

512 3 

1800 

32 

12.2 


We refer to spatial dissipative structures in a given 
time snapshot as states; they are identified as spatially- 
connected sets of points with current densities |j (cc) | ex¬ 
ceeding a fixed threshold, jthr 0. Typical states are 
thin, ribbon-like current sheets aligned with the z direc¬ 
tion. They occupy a small fraction of volume but ac¬ 
count for a large fraction of the overall resistive energy 
dissipation. Their lengths and widths span the inertial 
range, while thicknesses are localized inside the dissipa¬ 
tion range 13. l3|. 

Our present work extends this procedure into the tem¬ 
poral realm by applying a similar threshold criterion to 
the 4D space-time field j(x,t). We refer to the result¬ 
ing 4D spatiotemporal structures as processes or flare 
events. Our numerical algorithm first identifies the states 
in each snapshot as described above. It then connects the 
states in time by finding, for each given state, any other 
states in the two adjacent (past and future) snapshots 
having points connected to the same spatial region. In 
general, the states in one snapshot are not in bijective 
(i.e., one-to-one) correspondence with the states in the 
adjacent snapshot. This is due to interactions between 
structures, including mergers and divisions, and also the 
spontaneous formation of new states and destruction of 
old states. We refer to a sequence of states representing 
bijective evolution of a structure, beginning and ending 
with interactions, formation or destruction, as a path. 
Processes are then obtained as sets of paths connected 
via interactions. The natural and conservative approach 
for a temporal analysis is to study processes rather than 
individual paths, which become ambiguous upon inter¬ 
acting. 

We characterize processes by the following quantities. 
The process duration r is the time between the final 
state of the process and its initial state (normalized to 
Teddy)- The length L for a state is the maximum dis¬ 
tance between any two constituent points (normalized 
to the perpendicular box size L±). This is generalized 
for a process as the maximum length, L max , among con¬ 
stituent states. The instantaneous Ohmic energy dissi¬ 
pation rate for a state is £ = f dVrjj 2 , with integration 
over the constituent points of the state (normalized to 
average total energy dissipation rate, £tot ~ 1)- This 
is generalized to the total dissipated energy of a pro- 
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cess, E = J dt f dVr/j 2 , by combining the energy dis¬ 
sipation rates of all constituent states (normalized to 
E tot ~ £tot r eddy ~ 1). We also consider the peak energy 
dissipation rate, £ max , which is the maximum energy dis¬ 
sipation rate among constituent states. 


Results .— For the following analysis, we choose 
jthr/jrms ~ 6.8. This threshold is high enough to avoid 
percolation of processes through space and time. Pro¬ 
cesses that exist during the initial or final snapshots do 
not significantly contribute to the statistical results; we 
retain them in distributions for better statistics. 

An example process, with duration r rs 0.5 and 31 
distinct paths, is shown in Fig. [T] Representative states 
are shown (in green) on a subdomain of the simulation 
grid. We also show a schematic of the paths and interac¬ 
tions in the process. The process includes a division after 
the structure is stretched. A large number of paths are 
produced during the final stages, as the process decays 
toward the threshold. 

We now consider the statistical properties of the 
processes from the four simulations in Table HI The 
mean number of states per snapshot is (Wtate) = 
{194,288,657,1328}. For fixed cadence of A t~ l = 32, 
the mean number of processes per eddy turnover time 
is N pmc = {914,1271,4272,11608}, strongly increasing 
with Re. The most complex processes have ~ 10 3 con¬ 
stituent paths. We find a consistent asymmetry in the 
interactions: there are more divisions than mergers, with 
a ratio N mei /Ndi v = {0.84,0.78,0.80,0.82}. 

We show in Fig. [2] the probability distributions for 
dissipated energy E and for peak energy dissipation 
rate £ max . The distribution for dissipated energy, P(E), 
has a power-law tail with an index near —1.75±0.1, which 
is close to the analogous observations for total energy re¬ 
leased in solar flares 0, 0]- The power law extends 
across three orders of magnitude in E, from E rs 10 -5 
up to about E rs 10 -2 . For smaller E, the distribution 
is shallower and apparently non-universal, likely due to 
dissipation-range effects and threshold effects. With in¬ 
creasing Re , the power law extends to smaller E , consis¬ 
tent with the longer inertial range. The distribution for 
peak energy dissipation rate, P(£ max), exhibits a power 
law with index close to —2.0 ± 0.1 from f max « 10 -4 
to £max ~ 10 -2 . Similar indices are observed in distri¬ 
butions for peak hard X-ray flux in solar flares [e.g.,|44] 
and for energy dissipation rates £ of states [14]. 

The distribution for process durations r is also shown 
in Fig. [2] The durations extend to well above an eddy 
turnover time, sometimes comparable to the analyzed 
time interval. The distribution from t k 0.2 to r « 8 
can be fit to a power law with index near —3.2 ± 0.2, 
somewhat steeper than the indices ran ging between —2.2 
and —3.0 for solar flare durations 23l 32.144]. although 
close to the index —3.4 for rise times [281 ] . 

The process characteristics are related by strong cor¬ 


relations, with examples shown in the scatterplots in 
Fig. [3] We find that L ma x_~ t, previously inferred 
in solar flare observations [45[, while 5 max ~ r 2 and 
E ~ t 3 . These scalings are consistent with the estimate 
E = f dt f dVr)j 2 ~ r£ max ~ TV ma , x r)j 2 hl ~ t 3 , assum¬ 
ing that volume scales as length squared [H, [y[ and 
current densities are near }thr- From these correlations, 
E ~ £max, and hence P(E) is shallower than P(£ m ax ). 

Finally, we consider the evolution of processes over 
their durations, based on the constituent states at 
the given times. We focus on the energy dissipation 
rate, £{t) for 0 < t < r. Although this time-series 
is irregular for any given process, we find that pro¬ 
cesses of all durations exhibit similar average evolu¬ 
tions. In Fig. [4] we show the energy dissipation rate 
normalized to the peak, £(t )/£ max , versus normalized 
time, t/r, averaged for all processes. We find that 
(£ (t / t) / £ max ) ~ sin(7rt/r). Minor deviations suggest 
a temporal asymmetry, quantifiable by the first moment, 
{t/r) £ = f*(t/T)£(t)dt/ fo£(t)dt, which is 0.5 for sym¬ 
metric functions. We measure a very small but consis¬ 
tent asymmetry: for Re = {800,1250,1800}, (t/r)f = 
{0.483,0.483,0.476}. Hence a flare event grows slightly 
faster than it decays; this is similar to observed solar and 
stellar flares, although the asymmetry is amplified by the 
Neupert effect 28l 129. 32.?36. 46|. 


Conclusions .— In this Letter, we investigate the com¬ 
bined temporal and spatial intermittency of energy dissi¬ 
pation in numerical simulations of MHD turbulence. The 
conclusions are robust with respect to cadence, resolu¬ 
tion, and threshold. We find that a significant fraction 
of energy dissipation occurs in current sheets that are 
temporally organized in intense, long-lived flare events 
with durations that may span several large eddy turnover 
times. The process duration scales proportionally to its 
maximum length. The energy dissipated in these in¬ 
tense processes is distributed as a power law with in¬ 
dex near —1.75, implying the dominance of large, intense 
flare events. This can be compared to the spatial struc¬ 
tures at fixed times, which have a distribution of energy 
dissipation rates closer to the critical index of —2, sug¬ 
gesting that structures of all intensities instantaneously 
contribute equally to the overall energy dissipation rate 
T4]. 

We find that the distributions and scalings are in¬ 
sensitive to Re, suggesting universality. This is consis¬ 
tent with the fact that the structures are coherent across 
inertial-range scales, with only their thickness being set 
by the dissipation mechanism. The dependence of the 
statistics on Re has likely saturated at the relatively low 
Re considered here, making the results relevant for large 
Re in space and astrophysical turbulence. 

We find asymmetry in the divisions and mergers of 
current sheets, as well as in the evolution of the energy 
dissipation rate of a process. This temporal asymmetry 
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FIG. 1. States in a typical process with duration r rs 0.5, shown in green on a piece of the simulation lattice (with dimensions 
O.lOLx x 0.14Lx x 0.90Lx, without accounting for elongation of the lattice vertically). Also shown is a schematic of paths and 
interactions in the process, with red lines marking the times corresponding to the shown states. 



FIG. 2. The probability distributions for dissipated energy E, peak energy dissipation rate £ max , and duration r for Re = 800 
(red), Re = 1250 (blue), and Re = 1800 (green). 


can possibly be linked to the direct cascade of energy 
from large to small scales, giving large structures a ten¬ 
dency to divide into smaller structures. In future studies, 
temporal asymmetry may be a useful diagnostic for cur¬ 
rent sheet instabilities, including the tearing instability 

iffl- 


The present work lays out the foundation for a com¬ 
prehensive statistical analysis of dissipative processes in 
MHD turbulence, to appear in a future paper [5l|. The 
methodology can be applied to many other systems, in¬ 


cluding hydrodynamic turbulence M: line-tied MHD 

[HH |. avalanching 


53[ . kinetic plasma turbulence [54 
systems [Htij , and other complex dynamical systems. 
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FIG. 4. The evolution of energy dissipation rate relative to 
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durations, for Re = 800 (red), Re = 1250 (blue), and Re — 
1800 (green). This is well fit by a sine function (black), with 
small asymmetric deviations. 
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